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Abstract. We study the dynamics of vehicular traffic in a loop using a car-following model with the 
consideration of volume exclusions. In particular, we solve the steady state for the single-cluster case and 
derive fundamental diagrams, exhibiting two branches representative of entering and leaving the jam, 
respectively. By simulations we also observe that the speed average over all vehicles reaches the same 
constant at the steady states, regardless of the final clustering state. The autocorrelation functions for the 
overall speed average and single vehicle speed are investigated, each revealing a robust time scale. The 
effects of noises in vehicular acceleration are discussed. 

PACS. 89.40.Bb Land transportation - 05.40.Ca Noise 



1 Introduction 

The framework of statistical mechanics has reached ac- 
claimed triumph over the last century. The remarkable 
triumph simply lies on its success for reducing complicated 
many-body problems of Hamiltonian dynamics into simple 
thermodynamic relations. This success leads to specula- 
tions of applying its framework towards other many-body 
problems, such as the granular system [TJIMS], at which dis- 
sipation is introduced into the microscopically Newtonian 
dynamics. 

Meanwhile, a great deal of attention has been paid on 
self-driven sy stems @][5][B], with which it is much harder 
to build up the analogy to the framework of statistical 
mechanics. The traffic problem may provides an intrigu- 
ing example [511511 1 II 1 1 II 1 2 II 1 3 flUl 1 5 II 1 6 UTT] . In addition to 
the nonlinearity in its microscopic interactions, there are 
neither energy nor momentum conservations in a traffic 
system. Furthermore, the Newton's third law fails also, as 
each vehicle often makes responses according to the be- 
havior of its neighboring vehicle in front instead of the 
one behind. These features makes the problem almost in- 
accessible from the routes of traditional statistical me- 
chanics. Nevertheless, the studies and observations from 
traffic dynamics reveal some features commonly shared 
with statistical mechanics, such as phase transitions and 
fluid-dynamic features. 

In this work we aim to study the asymptotic dynam- 
ics and its dynamical fluctuations in a closed-boundary 
region from our traffic model. It is worth noting that ex- 
periments have been performed with such a circular road 
setting by real vehicles [TB]. The asymptotic evolution can 
be compared to a thermodynamic equilibrium in the lan- 
guage of statistical mechanics, as the overall pattern re- 
mains steady while the 'microscopic' objects keep on mo- 
tioning if possible. For the traffic system, the asymptotic 



regime bears time-translational symmetry apart from mi- 
nor fluctuations. Therefore macroscopic quantities such as 
the overall energy and momentum reach rather steady val- 
ues at this regime. 

The rest of this paper is organized as follows. In the 
next section we introduce the algorithms of our traffic 
model. The simulation results for the noise-free case are 
shown in Sec. [31 as it is shown that the system self-assembles 
into asymptotic steady states of various jamming patterns. 
Prompted by the observation of a cyclic, repetitive struc- 
ture in vehicular profiles, we solve for the steady state 
directly from the algorithms in Sec. [H and derive fun- 
damental diagrams and the retarded time thereby. The 
effects of noises in vehicular acceleration are studied from 
the autocorrelation functions introduced in Sec. [5[ and 
there turns out to be two robust timescales in the au- 
tocorrelation functions . Our results are summarized in 
Sec. EJ 



2 Model 

In our current work, we choose to use a continuous car- 
following model instead of a cellular automaton mo del [XT] . 
because the former gives a better description about dy- 
namics. In the simplest car-following model, the propor- 
tionate control [5], the acceleration of a vehicle is depen- 
dent on its velocity difference to the leading vehicle. How- 
ever, based on our daily driving experiences, this corre- 
lation decays as its headway (the distance to the leading 
vehicle) increases, and at a large distance the vehicle tends 
to accelerate towards some optimal velocity vq. Therefore 
we adapt a modified car-following model that interpolates 
between these two scenarios. Moreover, we put in an op- 
tional noise term for the acceleration of each vehicle, rep- 
resenting internal fluctuations of vehicles. 
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As to our algorithm, the acceleration of the ith vehicle 
in free traffic is: 



where 



v i+1 + (v -v i+1 )(l-e- Ax */ D f): 



(1) 



(2) 



as Df is the characteristic car- following distance, Axi is 
the distance to the leading vehicle, and r\ is an optional 
noise term representing the internal fluctuations of each 
driver. In addition to the modified car-following algorithm 
as Eq. [TJ we put on an extra braking algorithm to enforce 
the volume-exclusion criterion Axi > D c , where D c is the 
length of each vehicle, as in this work we assume the ve- 
hicles to have a uniform size. Thus if at one moment the 
profile shows Axi — D c for some i and the projected up- 
date from Eq. [TJ results in a violation of volume exclusion, 
then we reset Vi — 0. This algorithm mimics our driving 
experience in the sense that a driver tends to stop to avoid 
an accident if the leading vehicle is in very close proxim- 
ity. At last, we require that a vehicle at stop will start to 
accelerate according to the algorithm in Eq. Q] only if its 
distance to the leading vehicle is larger than some safety 
distance D s . 

From the accelerating algorithm in Eq. [TJ our model 
bears resemblance to the intelligent driver model [13], which 
also takes into consideration the dependence on both the 
velocity difference and distance to the leading vehicle. 
Furthermore, the consideration of a noise term and the 
volume-exclusion criterion is characteristic in the Nagel- 
Schreckenberg model [TTJ. Thus one can treat our model 
as a continuous version of the Nagel-Schreckenberg model 
with more emphasis on the acceleration and deceleration 
dynamics. Also the enforcement of the volume-exclusion 
criterion points out an extra decelerating mechanism which 
is much faster than the usual car-following mechanism, al- 
though in our model this extra deceleration is performed 
instantly. 

In this work we apply a periodic boundary condition, 
so that the vehicles can be treated as proceeding on a ring 
of length L. Unless mentioned we set the vehicle length 
D c =3 meters, safety distance D S =Q meters, the character- 
istic following distance D/=60 meters, and £=1000 me- 
ters. We set the optimal velocity vq = 25m/s = 90km/hr. 
The parameter A represents the pace of car-following adap- 
tation, and the larger A is, the sooner the vehicle i adjusts 
to vf cxt . In this work we set A = 0.15 s _1 , such that it 
takes about 20 seconds for a vehicle at stop to reach ve- 
locity 0.95vo if its front is empty of other vehicles. During 
our simulations we set the time interval between successive 
updates to be At = 0.001s. We use the SI units through- 
out this report. As to the initial conditions for all of our 
simulations, we set the vehicles to distribute uniformly 
over the road, while we use various speed profiles. 
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Fig. 1. Time evolution for n st0 p , the total number of stopped 
vehicles, and « ave , the speed average over all vehicles. The 
curves show results of trials over different initial speed dis- 
tribution with the setting D s = 6m. The curves are labeled by 
the number of jammed clusters observed at their asymptotic 
steady states. 



3 Simulation results: noise-free case 

Fig.rfjshows the time evolution in our simulations for var- 
ious initial conditions, as n st0 p is the number of vehicles 
stopped and v avc is the speed average over all vehicles. Al- 
though it is well known that the algorithm of proportion- 
ate control does lead to stable traffic asymptotically [9], 
the algorithm itself does not guarantee that vehicles do 
not overtake their front neighbors. With the extra con- 
sideration of volume exclusion, jams occur due to the 
propagation of the sudden-brake behavior. For our closed- 
boundary system, the propagation of jamming often comes 
to a stop upon reaching a depleted zone, as the traffic dy- 
namics finally evolves towards an asymptotic steady state 
with some rather stabilized pattern. At the asymptotic 
steady state the quantities n stop and v avc exhibit slight, 
regular fluctuations only. In our simulations we have not 
observed any chaotic speed profiles asymptotically for the 
noise-free case. Fig. [5] shows an example in which the 
asymptotic steady state consists of five jammed clusters. 
One observes that at large time the vehicles tend to re- 
peat a single pattern cyclically, as each vehicle repeats the 
behavior of its leading vehicle with a fixed time delay. 
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Fig. 2. The trajectories of vehicles during one simulation. At 
large time the vehicles tend to repeat a single pattern cyclically, 
as each vehicle repeats the behavior of its leading vehicle with 
a time delay. 



It is worth noting that Fig. Q] shows the speed average 
w a ve reaches approximately the same asymptotic value de- 
spite the various initial conditions considered. Therefore 
no matter how many jammed clusters it forms, the speed 
average of the system stays the same asymptotically, apart 
from minor fluctuations. The only exception happens for 
the non-jammed state, in which w avc = vq. This is less in- 
tuitive as we often think that v aV c would be smaller if more 
jams are encountered at the traffic. This fact can be helped 
understood partially noting that in Fig. [TJ one finds less 
total number of stopped vehicles instead as there are more 
jammed clusters. However, it remains to account for the 
remarkable coincidence that D aV e reaches the asymptotic 
same value regardless of jamming patterns. 

As a comparison, we set the safety distance D s to zero 
and find that D avo does not reach a fixed value over dif- 
ferent jamming pattern formations. Fig. [3] shows the time 
evolution of n st0 p and i> avo for two examples. Furthermore, 
we also find that it takes a much longer time for the fluctu- 
ations to die out before the system reaches an asymptotic 
steady state. Also the lack of a safety distance implies that 
large fluctuations could arise at a slight perturbation. For 
example, if we set that initially Vi(t = 0) = vq for all i 
except vi = vq — 20m/s. While for the case with D s = 6m 
we observe a single cluster pattern, a four-cluster pattern 
emerges for the case D s = 0m, which suggests that with- 
out a safety distance, subsequent jams can be triggered 
easily from a local disturbance. 



4 Steady-state solution of single-cluster 
formation 

From our simulation results for the noise-free case, we 
find at as time proceeds, the traffic evolves towards a 
steady pattern, as is evidenced by the stabilized macro- 
scopic quantities, such as the averaged speed and number 
of stopped vehicles. Furthermore, the position profile of 
vehicles reaches a cyclic pattern asymptotically. Assuming 
that asymptotically all the fluctuations and irregularities 
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Fig. 3. Time evolution for n stop and u ave for the case D s = 0m. 
The curves show results of trials over different initial speed 
distribution, as curves are labeled by the number of jammed 
clusters observed at their asymptotic steady states. 



turn to die out, and vehicles are proceeding with same 
velocity pattern except with time shifts, one can solve for 
the asymptotic behavior from our car-following algorithms 
described in Sec. [2j 

Let us define t — i m ; n to be the time that the ith 
vehicle starts to dissolve from the jam, and t — to be 
the time for it to reenter the jam. Therefore one has 



Ax(t = 0) = x i+1 (0) - Xi(0) = D c , 



and 



Ax(t = t n 



(3) 



(4) 



From the cyclic pattern we assume that at the steady state 
the velocity profile of each vehicle shifted from that of its 
leading vehicle simply by a time delay r, i.e., 



«»(*) = v i+ i(t - t) . 
By putting Eqs. [2] and [3] into Eq. [1] we get 

1 dvj 



(5) 



A dt 



= v l+1 (t) - Vi (t) + [v - v i+1 (t)} ■ [1 - e-^W/^] 



= Vq - Vi(t) - [vq - Vi(t + t)] 



exp 



Vi(t')) 



D 



f 



,(6) 
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where we have used the fact that 



Axi(t) 



dt' [vi+i(lf)-Vi(!?)] + Axi(p) 



= / dt'[Vi(t' + T)-Vi(t')]+D c 





t + T rT 

dt'vi(t') - / dt'vi(t') + D c 
Jo 

t + T 

dt' Vi(t')+D c . 

t 



(7) 



Note that we assume that Vi(t) = for t > 0. This as- 
sumption holds as long as each vehicle comes to a stop 
before its leading one starts. Eq.|5]is an integro-differential 
equation that cannot be solved by simple numerical inte- 
gration. Since the delayed time r is itself yet to be known, 
we resort to an iteration procedure instead. To determine 
the value of r, we need an extra identity, which is given 
by Eq. [4] (along with the use of Eq. EJ: 



AXi (tmiii ) D s 



dt'vi(t') + D c 



(8) 



Note that the boundary condition for Eq. [6] is Vi(t = 
tmin) — 0. The value of vi before reentering the jam Vi(0~) 
will be determined upon solving Eq. [5] Thus each itera- 
tion consists of two steps. First we decide a new value of 
t to match the identity in Eq. [5] using the previous speed 
profile. Then we derive the updated speed profile accord- 
ing to Eq. [5] In practice we find it very effective to rewrite 
Eq.i] as 

e- At ^[e At fa(t)-«o)] = \[vi(t+T)-vo]e-A**W D t , (9) 
at 

and for each iteration we plug in the old speed profile on 
the right-hand side and derive the updated profile from 
the left-hand side of Eq. O By this algorithm we can get 
a steady-state solution with the sums of residual squares 
less than xlCP 4 within 15 iterations, as we set the time 
interval between successive updates to be At — 0.001. 
Before iteration starts we set our initial guess of the speed 
profile to be Vi (t) = v -{l— cxp[— X(t— t min )]}-[l~ exp(Ai)]. 

In Fig. [4] the steady-state solutions of the speed pro- 
files for the cases i m j n = —40, -20, and -10 seconds. We 
stress again here that they are derived from the single- 
cluster steady-state condition. Note that — i m i n indicates 
the time of free traffic for a single vehicle between succes- 
sive stops. One can observe a sharp, discontinuous change 
in the acceleration of the vehicle, which occurs at the mo- 
ment when its leading vehicle stops. Furthermore, there 
exist great asymmetry between the accelerating and the 
decelerating parts of the speed profile, a result of the fea- 
ture that there is no time-reversal symmetry in the micro- 
scopic traffic dynamics. The second discontinuity in Fig.H] 
happens at t = 0, when the speed drops to zero abruptly. 
This implies that the car-following algorithm we used in 
Eq. [T] cannot guarantee a collision-free environment and 
the driver has to resort to an extra, instantaneous braking 
behavior to maintain the volume exclusion principle. 




Fig. 4. Speed curves of steady state traffic for i m i n = —40, 
-20, and -10 seconds. 
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Fig. 5. Fundamental diagrams of steady state traffic for t mlu 
-40, -20, and -10 seconds. 



From the steady-state solution we can obtain funda- 
mental traffic diagrams as follows. First of all one has 
j = pv, where j is the traffic flow and p is the density of 
vehicles. Since we already have the speed profile Vi(t), we 
can use the definition that p(t) — 1/Axi(t), the inverse of 
the spacing between successive vehicles. Therefore we get 
j(t) = Vi(t)/Axi(t), as the fundamental diagram is plot- 
ted in Fig.[5]by collecting the serial data of (p(t),j(t)). For 
each t m i n the fundamental diagram shows two branches. 
The upper branch represents the dynamics towards a jam, 
while the lower one shows the dynamics coming out of 
the jam. The results shows that the slope for the upper 
branch in the fundamental diagram changes significantly 
at small |t m i n as t m - m modifies, while at larger |t m i n | the 
curves merge toward a single one. Note that since the lower 
branch represents dynamics coming out of a jam, where 
the vehicles are accelerating greatly from stop, the curves 
at this branch are more sensitive to the definitions of j 
and p. 

One can estimate the receding velocity of the jamming 
front Wj am from the steady-state solution via two routes. 
First, realizing that for each time interval r, a vehicle en- 
ters into the jam while another leaves the jam, one con- 
cludes that Uj am = —D c /t. The second method follows 
from the derivation of the shock- wave velocity in Ref. [S] . 
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Hence Vj am can be estimated from Fig. [S] by the slope of 
the line connecting the points of the jammed state (where 
j = 0) and the free state (estimated by the point with 
lowest density allowed in the fundamental diagram). For 



the case t„ 



-40 the both methods give Vj, 



1.11 



-^free — 



dtVi(t) - Vj a 



(10) 



due to the fact that the jamming front moves backwards 
with velocity i>j am while the vehicle proceeds. Alterna- 
tively, one can derive Lf ree by summing over the spacings 
between successive vehicles at t = t min : 

n n 
Lfr cc = ^ (*min) = ^ ^ min + JT ) ' 



where n is the number of vehicles at free traffic: 



"-fro 



T 



(12) 



Thus for each t m j n one gets a corresponding Lf ree and 
Ufree at the single-cluster steady state. Through this one- 
to one correspondence we can treat £f 1G o as a function 
of nf r ee- Therefore if n vehicles are running on a loop 
of circumference L > Lf ree (n), then one cannot observe 
any single-cluster formation at the steady state. To verify 
this argument, we first note that the free-traffic section 
in single-cluster formations bears the property that the 
residual (empty) space at free traffic is a monotonically 
increasing of n: 

£free(«) - nD c > L {ree (n') - n D c if n > n' . (13) 

Now assume n vehicles are distributed at a loop with cir- 
cumference L > Lf ree (n). If at the steady state it forms a 
single cluster with n — n' jammed vehicles, then the length 
of the free-traffic section is L — (n — n')D c . Therefore we 
have 

L {iee (n') = L-(n-n')D c . (14) 
However, Eq. [14] implies that 

Lhcc{n') - n'D c = L - nD c > L trcc (n) - nD c (15) 

for some n' < n, which violates Eq. 1131 Hence we con- 
clude that for n vehicles distributed at L > Lf vcc (n) single- 
cluster formation does not exist at the steady state. Fig. [5] 
shows such a phase diagram while Lfr eo is evaluated us- 
ing the second method. The results from the first method 
follow very closely, with residuals less than the size of a ve- 
hicle. Note that the condition that L < L[ YCC (n) does not 
necessarily imply the existence of a jammed steady state, 
as states with a homogeneous speed distribution Vi = vq 
also act as stable solutions. 
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m/s. Furthermore, the value Wj am changes less significantly 'g 
over t m i n . For example, our solution shows Wj am ~ — 1.10 j 
m/s when t m - ln = —10s. 

The length of the free-traffic section at the steady 
state Lfree can be evaluated in the following manner. Since 
—t m in is the duration of time for each vehicle between suc- 
cessive stops, one has 
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Fig. 6. Phase diagram for single-cluster formation. L is the 
circumference of the loop, and n is the total number of vehicles. 



5 The noise-added case: 
functions 



autocorrelation 



Instead of putting random noises on all drivers at all times, 
we note that noises occur rather occasionally and choose 
put the noise term in our model as follows. During each 
time step At, each vehicle has a probability p to suffer from 
random noise rj in its acceleration due to disturbance or 
else. The random variable r] bears a uniform distribution 
ranging between — r/o and t/q. Throughout the course of 
simulations we use the value rjo = 1000m/s 2 . Although 
the value of rj looks rather unphysical, within the time 
between successive updates At = 0.001s we find that each 
noisy acceleration results in a change of a speed up to 
rjo At = lm/s, or equivalently, 3.6km/hr. 

For the noise-added case, the vehicles do not evolve 
into a cyclic, steady-state pattern asymptotically due to 
persistent internal fluctuations. Nevertheless, asymptot- 
ically the system still evolves into a stage resembling a 
thermal equilibrium. Hereby we study the nature of the 
asymptotic kinetics via autocorrelation functions. We de- 
fine the autocorrelation function for the speed average 

Cave(^) as 



Cave if) 



_ ((faveft') ~ (^avc))(favc(^ + t) - (^avc))) 
(«avo ~ (favo) 2 } 



and the autocorrelation for the speed of a certain vehicle 
d(t) as 



Cx{t) 



{Vi))(ViW+t)-{Vi))) 



(17) 



as the bracket represents time average. 

Fig- shows the autocorrelation function C ave (i) over 
various values of noise probability p. We find that for 
small p the relaxation shows short-time oscillatory behav- 
ior, with the oscillation period close to 2.7s. This oscil- 
lation period is approximately equal to the delayed time 
t, which we have derived from our steady-state analysis 
for single-cluster formation. Hence the short-time oscilla- 
tions represent the synchronized behavior of vehicles en- 
tering and leaving jamming sections. We find that for the 
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Fig. 7. Autocorrelation function C a ve(i) of speed average « a ve- 
The parameter p represents the likelihood of noise occurrence. 



noise-free case the same oscillation period persists as well 
for steady states of multiple clusters. Furthermore, when 
the probability of noise occurrence p increases, the oscilla- 
tory period also increases slightly. As p further increases, 
the synchronized behavior eventually diminishes, and the 
short-time oscillations become absent in the relaxation of 
w a vc, and the long-term memory possessed by the system 
in the noise- free case eventually gets destroyed, as C avc (t) 
reaches a short-time exponential decay at large p. 

The relaxation of single- vehicle speed fluctuation C\ (t) 
is presented in Fig. [51 For the noise-free case we observe 
peaks of height one for every time interval of approxi- 
mately 162s. By checking our single-cluster simulation re- 
sult we find this is the time required for a vehicle to repeat 
its speed profile at the steady state. Furthermore, by ex- 
amining our results over several noise-free cases we find 
the same time interval persists even for systems with mul- 
tiple jammed clusters. Meanwhile, for multi-cluster cases 
the autocorrelation function C\(t) exhibits smaller peaks 
in between. When noises are added, we find this corre- 
lation to be rather robust, compared with the collective 
relaxation of v ave . The peaks remain significant even at 
p = 0.01, while a more detailed check suggests the time 
interval between successive peaks slightly widens as p in- 
creases. 

For the noise-free case, the speed average v avc exhibits 
small, regular fluctuations due to the repetitive stopping 
and accelerating motions of vehicles at jamming bound- 
aries. When the noise level increases this regular fluctu- 
ation will be easily overpowered because of its small am- 
plitude. Meanwhile, the structure or the pattern of the 
traffic are less susceptible to the noises. It thus explains 
why the short-time oscillations in C ave (i) diminish quickly 
with noise increasing, while the major peaks in C\(t) re- 
main up to a relatively high noise level. However, it re- 
mains to understand why the time intervals between peaks 
stay unchanged or are only slightly modified over different 
jamming patterns and noise levels. 

One may note from Fig. [8] that for the noise-free case 
the autocorrelation functions C\ (t) exhibit symmetric pat- 
terns between successive largest peaks. This symmetry 
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Fig. 8. Autocorrelation function Ci(t) of single-vehicle speed 
Vi. The parameter p represents the likelihood of noise occur- 
rence. Upper: noise-free cases (p = 0); Lower: with added 
noises. The noise-free curves are labeled by the number of clus- 
ters formed at the steady states. 



results from the fact that at the asymptotic state time- 
reversal symmetry holds for the autocorrelation function: 
Ci(t) = Ci(— t), due to the fact that time-translational 
symmetry holds asymptotically: (v(t')v(t' + <)) = {v(t' — 
t)v{t')). And for the noise-free case the system reaches a 
steady state asymptotically with a cyclic pattern. There- 
fore d(t) = d{t + T) = C\(-t) =d(T- t), where T is 
the cyclic period. On the other hand, as the noises turn 
on, the autocorrelation function C\(t) fails to be strictly 
periodic, as asymmetric pattern emerges therein. 



6 Summary 

It is important to understand the many-body dynamics 
of self-driven objects, which appears so often in biologi- 
cal and even social systems. The highway traffic dynam- 
ics provides a simple platform towards understanding the 
possible common behavior of self-driven systems. Although 
the microscopic dynamic equations are highly nonlinear 
compared with Newtonian dynamics, the traffic system 
with closed boundaries evolves towards a steady state for 
the noise-free case. We observe from simulations that the 
speed average over all vehicles eventually reaches approx- 
imately the same value, regardless of the clustering pat- 
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terns. This remarkable feature is absent if we do not allow 
any safety distance before stopped vehicles start moving. 

The steady states can be solved knowing that succes- 
sive vehicles are following same speed profile with some 
retarded time r. And therefore one can obtain fundamen- 
tal diagrams for the noise-free case. Moreover, we have 
derived a criterion for possible single-cluster jam forma- 
tion. 

Meanwhile, as noises are introduced for individual ve- 
hicles, the asymptotic evolution of the system can be ana- 
lyzed via autocorrelation functions. The autocorrelation of 
the speed average C avc (t) exhibits rapid oscillations with 
a period characteristic of the retarded time r. Again this 
period is robust over different traffic patterns, and this 
short-time regular oscillations for the speed average will 
be overpowered as the likelihood for noise occurrence p in- 
creases. Furthermore, the single-vehicle speed autocorre- 
lation function C\ (t) shows periodic behavior with slower 
oscillations for noise-free cases. The period observed cor- 
responds to the time for a single vehicle to repeat its speed 
profile. Again this timescale is robust over our various sim- 
ulation trials. When noises are added we find these peri- 
odic peaks can live up to a high noise level. 



References 

1. S. F. Edwards and R. B. S. Oakeshott, Physica A 157, 1080 
(1989). 

2. H. J. Herrmann, Physica A 313, 188 (2002). 

3. Y. M. Jiang and M. Liu, Gran. Matt. 11, 139 (2009). 

4. T. Viscek, A. Czirok, E. Ben-Jacob, I. Cohen, and O. 
Shochet, Phys. Rev. Lett. 75, 1226 (1995). 

5. G. M. Schiiz, in Phase Transitions and Critical Phenomena, 
eds. C. Domb and J. L. Lebowitz (Academic Press, London, 
2000), Vol. 19. 

6. F. Schweitzer, Brownian Agents and Active Particles 
(Springer, Berlin, 2001). 

7. D. Helbing, Rev. Mod. Phys. 73, 1067 (2001). 

8. M. J. Lighthill and G. B. Whitham, Proc. R. Soc. London, 
Ser. A 229, 317 (1955). 

9. R. E. Chandler, R. Herman, and E. W. Montroll, Oper. Res. 
6, 165 (1958). 

10. I. Progigine and R. Herman, Kinematic Theory of Vehic- 
ular Traffic (Elsevier, New York, 1971). 

11. K. Nagel and M. Schreckenberg, J. Phys. I France 2, 2221 
(1992). 

12. M. Bando, K. Hasebe, A. Nakayama, A. Shibata, and Y. 
Sugiyama, Phys. Rev. E 51, 1035 (1995). 

13. M. Treiber, A. Hennecke, and D. Helbing, Phys. Rev. E 
62, 1805 (2000). 

14. B. S. Kerner, The Physics of Traffic (Springer, Berlin, 
2004). 

15. R. Mahnke, J. Kaupuzs, and I. Lubashevsky, Phys. Rep. 
408, 1 (2005). 

16. Y. Sugiyama et ai, New J. Phys. 10, 033001 (2008). 

17. B. S. Kerner, Introduction to Modern Traffic Flow Theory 
and Control: The Long Road to Three-Phase Traffic Theory 
(Springer, Berlin, 2009). 



